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Evidence for a hot dust-free inner disk around 51 Oph ^ 
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Abstract. We report on the observation of CO bandhead emission around 5 1 Oph (Av = 2). A high resolving power {R ^ 10,000) 
spectrum was obtained with the infrared spectrometer ISAAC mounted on VLT -ANTU. Modeling of the profile suggests that 
the hot (Tjas = 2000-4000 K) and dense (riu > 10'" cm"^^) molecular material as probed by the CO bandhead is located in 
the inner AU of a Keplerian disk viewed almost edge-on. Combined with the observation of cooler gas (rgas = 500-900 K) 
by ISO-SWS and the lack of cold material, our data suggest that the disk around 5 1 Oph is essentially warm and small. We 
demonstrate the presence of a dust-free inner disk that extents from the inner truncation radius until the dust sublimation radius. 
The disk around 51 Oph may be in a rare transition state toward a small debris disk object. 
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1. Introduction 

Disks are commonly found around low- and intermediate-mass 
■ pre-main-sequence stars and are understood as natural byprod- 
ucts of the conservation of angular momentum during star for- 
mation. The structure of the inner few Astronomical Units 
(AU) of protoplanetary disks plays an important role in the 
transfer of matter onto the star and the formation of terres- 
trial planets. The temperature in the inner AU of disks can at- 
' tain a few thousand Kelvin and the density can be greater than 
lO'^ cm"^. Those conditions are required for CO bandheads 
(Av = 2) to emit. Indeed, CO bandhead emission h as been de- 
tected aro und young low- and high-mass stars (e.g.. lBik & Thil 
I2004 H and'N aiita et aljEoOOl for an overview). The inner disk 
is likely gas-rich but dust-free since dust grains evaporate at 
T ^ 1500 K. 

51 Oph is a peculiar Be (B9.5IIIe) star located at 130 pc 
with an observed ro tation velocity vsin/ - 267 ± 5 kms"' 
jPunkin et alJll997l) : however, the exact spectral type is still 
disputed. The Spectral Energy Distribution (SED) shows a 
strong near- till mid- infrared excess but contrary to its lower- 
mass counterparts the Herbig Ae stars, the flux in the (sub)- 
millimeter domain is low, indicating a lack of small cold 
dust grains. 51 Oph is the only Be star in a sample of 101 
whose infrared excess is dominated bv w arm dust and not 
by free-free emission dWaters et alJ Il988{) . Both fits to the 
SED and direct imaging of the disk around 51 Oph agree 
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* Based on observations collected at the European Southern 
Observatory at La Silla and Paranal, Chile (ESO Programme 68-C- 
0474) 



on the compactness of the disk Jjavawardhana et alJ l200lt 
iLeinert et alJ2004tlWaters et al.ll988h . The disk probablv does 
not extent beyond 100 AU. Silicate emission features were 
first detect ed by Faiardo -Acosta et al. (1993) and further an- 
alyzed by iBouwman et alJ (1200 ll) . The UV and optical part 
of the SED are well fitted by a Kurucz model and show no 
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Ima). However, in- 



is low 

falling ionic and atomic gases were observed which prompted 
Gradv et al. ( 1991) to claim th at 5 1 Oph is in the same evolu- 
tionary state as y6 Pic (see also lRoberge et al]l2002h . Contrary 
to y6 Pic, large amount of molecular gas (CO, CO2, H2O and 
NO) has been seen around 5 1 Oph by the Short- Wavelength- 
Spectrometer on board the I nfrared Space Observatory (ISO) 
(I van den Ancker et al.l l200lh . A controversy remains about 
the quantity of CO. Absorption lines by C O in the Far - 
Ultraviolet were not detected by FUSE jRoberge et alJl2002h . 
Nevertheless, all studies suggest that the circumstellar matter 
is most likely in the form of a Keplerian disk rather than in 
a spherical shell. Finally, the fractional excess IR luminosi- 
ties Lir/L, of 51 Oph is 0.028, a value between the younger 
Herbig Ae stars (>0.1) and the Vega-like objects (10^^-10"^). 
We obtained a high resolution A'-band spectrum of 5 1 Oph cen- 
tered around the position of a CO bandhead emission in or- 
der to clarify the detection of large amount of CO (N{CO) - 
(0.1 - 10) X 10^' cm"^) by ISO and to characterize the struc- 
ture of the inner dust-poor molecular disk. After describing the 
observation and data reduction in Section|2l a fit by a synthetic 
spectrum emitted by a Keplerian disk is shown in Section |3l 
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Fig. 1. Continuum subtracted CO bandheads toward 51 Oph in 
arbitrary units. Four bandheads are detected. The spectrum of 
the last two bandheads is heavily blended. 

The discussion in Sectionl^focuses on the structure of the in- 
ner disk around 5 1 Oph. 

2. Observation and Data Reduction 

The ^-band spectrum was obtained with the near-infrared fa- 
cility ISAAC mounted on the Very Large Telescope (VLT- 
ANTU) on March 15* 2002 at resolving power of 10,000 (slit 
width of 0.3") in Service mode. The spectrum toward 51 Oph 
is part of a survey of CO bandhead in Herbig Ae stars. The data 
were reduced in a standard way using IRAF. We made use of 
flatfield and arc frames taken by the ESO staff during the day. 
Standard stars of spectral type A observed at similar airmass 
than 5 1 Oph were used to correct for telluric OH emission and 
absorption. 
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Fig. 2. Fits to the CO first bandhead (v - 2 - 0). In the up- 
per panel the thin line corresponds to the normalized observed 
spectrum while the thick line is the best fit if the disk is viewed 
at an inclination of 45 ° . The lower plot is the residual of the fit. 
The lower panel shows the best fit obtained with a disk viewed 
at an inclination of 88 °. The complete set of parameters is re- 
ported in Table [2 



3. Results and l\/lodeling 

The continuum subtracted and normalized spectrum is shown 
in Fig.^ Four bandheads are detected with high signal-to-noise 
ratio (S/N>100), namely Av = 2 - (A - 2.2935 jum), Av = 
3 - I (A ^ 2.3227 fim), Av ^ 4 - 2 (A ^ 2.3535 fim), Av = 
5 - 3 (/I ^ 2.3829 /im). The first bandhead shows a prominent 
blue wing. The higher bandheads are contaminated by emission 
from the P-branch from lower bandheads. 

We decided to fit the first bandhead only since the emis- 
sion from higher bandheads is blended and contributions from 
different bandheads are difficult to disentangle. A spherically 
symmetric and radially expanding wind results in flat-topped 
profiles for optically thin lines (.Beals, .1931) . A flat-topped 
profile cannot reproduce the blue-wing seen in our spectrum; 
therefore our data definitely rule out hot gas in an expand- 
ing shell. The signature of CO bandhead emission from a nar- 
row Keplerian rotating disk is commonly found around young 
stars over a large range of masses. Therefore, we constructed 
a simple model that is detailed below. The distribution of CO 
lines generated by a plane-parallel isothermal slab is convolved 



with the double-peaked profile produced by a naiTow ring ex- 
tending from 7?n,in to /?niax- The ring rotates at Keplerian ve- 
locity v'Kep - VCMT/t. T he line frequencies an d Einstein 
coefficients are provided by IChandra et alJ ( Il996l) . The rota- 
tional levels within each vibrational level are assumed to be 
in local thermodynamic equilibrium because the critical den- 
sity for thermalization of the rotational levels is relatively low 
(nrot,ciit - 10^ cm"^). On the other hand, the vibrational lev- 
els are only thermalized at densities greater than 10'" cm"^. 
Therefore, each vibrational level can be characterized by its 
own vibrational temperature Tyn,. A large range of tempera- 
tures was explored (T- 1000-4000 K) and the maximum tem- 
perature allowed for the dust is set at 1500 K. The gas tem- 
perature was allowed to vary with a power-law of index - 
3/4 (T(r) = Tni r/Rminy^^'^), typical for a flat accretion disk 
(lAdams et alll988i) and the gas surface density decreases with 
radius as r"' (N(r) = Afo('"/^min) ') ■ The free parameters of the 
models are the inclination /, the inner and outer radius (/?min and 
^max) of the hot CO emitting area, and the temperature and CO 
column density at the inner radius (To and A^o). A grid of mod- 
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els was performed and the model with the lowest reduced chi- 
square is shown in the lower panel of Fig.|3 The gas is located 
in a Keplerian disk seen with inclination angle of 88^35° with 
respect to the rotational axis. The CO column density and tem- 
perature at 0.15 AU are (0.17-2.5) lO^" cm^^ and 2850±500K 
respectively. The amount of CO gas can be inferred from the 
modeling and is M(CO) - (0.17 - 2.5) 10"'° Mq. Our results 
show that the CO gas is located well within the inner first AU 
from the central star (Table The maximum distance R^ax 
where the CO bandhead is emitted depends on the inclina- 
tion, which is not a well constrained parameter, but is probably 
smaller than the dust sublimation radius because the gas is at 
a much higher temperature than the dust sublimation tempera- 
ture. 

Table 1. Best fit parameters. / is the inclination, and 7?,nax 
are the emission area inner and outer radius, ro(CO) is the gas 
temperature at and No{CO) is the column density of CO at 
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-^min 

(AU) 


(AU) 


T„(CO) 
(K) 


A'o(CO) 
(cm-2) 




0.15±0.05 


0.35±0.05 


2850±500 


(0.17-2.5) 10^° 



4. Discussion 

The observation of CO bandhead confirms the large column 
density of CO seen by ISO-SWS. Moreover, the profile of the 
emission lines indicates that the gas is seen almost edge-on. 
If the disk around 51 Oph is seen with an inclination of 88^35 
and the disk is geometrical thin (flat disk), then the ISO-SWS 
and the FUSE data can be reconciled. In the FUSE absorption 
study, the line-of-sight does not intercept the major part of the 
disk while emission studies with ISO or ISAAC are not sen- 
sitive to the actual projection angle /. On the other hand, hot 
gas of N I, S II and Fe iii can be detected in absorption against 
the stellar photospheric emission because the evaporating bod- 
ies may follow a trajectory that is inclined with respect to the 
equatorial plane jBeust et alj200ll) . 

In the absence of extinction by dust grains, the principal de- 
struction agent of CO is photodissociation by stellar ultraviolet 
photons. However, the CO column densities (10^*^ - 10^' cm^^) 
are well above the required va lue for the CO mo lecules to 
self-shield (A^(CO) ~ 10'^ cm'^ lvan Dishoeck & Black.1988.) . 
Moreover, at high-temperature and density the rapid forma- 
tion of CO molecules by the neutral-neutral reaction C-i-OH 
^CO-i-H can easily compensate for the photodestruction. The 
large abundance s in H7O, CO? and CO tes tify of a hot and 
dense chemistry ( Van den Ancker et alJ200ll) . 

It is interesting to compare the inner (i^min) and outer (/?,nax) 
CO emitting radius to other relevant radii in the disk. The 
co-rotation radius r^t - GM^/nl is defined as the distance 
from the central star where the st ar and the Kepler ian disk ro- 
tate at the same speed a (e.g.. lShuetaTlll994 . Assuming 



that Q, = vsin/ - 267 km s"' (i.e. / = 90°), we ob- 
tain Trot - 0.05 AU (Table |2j- Disk material falls on to the 
star surface only when the disk truncation, whose value is 
close to the magnetospheric radius r^, is smaller than the co- 
rotation radius for an accretion disk. The magnetospheric ra- 
dius depends on the stellar magnetic field and the disk accre- 
tion rate, which are unknown for 51 Oph. We therefore as- 
sume that rm - Ti-ot. Dust grains do not exist above their subli- 
mation temperature, which is around 1500 K. The dust sub- 
limation radius is the distance from the star beyond which 
dust can condense: = y/QuiL, + L^^c)l^6ncrlT^^^ where 
Qr - Qabsia, r»)/gabs(fl, Tsub) the ratio of the dust absorp- 
tion efliciency at stellar temperature T, to its emission effi- 
ciency at the dust sublimation temperature Tjub, a is the mean 
grain radius and cr is the Stefan constant. For large silicate 
grains (a > 1 fim), Qr is relatively insensitive to the stel- 
lar effective temperature and close to unity because most stel- 
lar radiation lies at wavelengths shorter than the grain radius. 
For smaller grains, the value of Qr is significantly increased 
dMonnier & Millan-G abet 2002*). We found for 51 Oph that 
rd = 0.56 AU (see Table |2} assuming Lace << - 260 
and Qr - I (large grains). If the grains were as small as 0. 1 ^va 
in radius, then Qr ^ 20 and ra = 4.50 AU. 

The disk around 5 1 Oph is not strongly accreting as testi- 
fied by the absence of UV excess. The disk is passively heated 
with the gas in the upper layer being warmer than the mid- 
plane. It should be noticed that the dust sublimation radius 
is a lower limit where dust grains can exist. Radiative pres- 
sure, which is particularly effective for a B9.5V star, will push 
the grains much further out than r^. At radii greater than the 
dust sublimation radius, the dust and gas are probably ther- 
mally coupled and emission in the near-infrared will be dom- 
inated by the dust, reducing drastically the gas emission lines 
over continuum contrast. The large difference between and 
Trot (0.5-2.45 AU) implies the presence of a gas-rich dust- 
poor inner disk extending from the truncation radius, which 
is close to the magnetospheric radius, till the dust sublimation 
radius. It should be noticed that beyond 1 AU the gas will be- 
come too cold to excite the high vibrational levels (v > 2). 
On the other hand high-7 CO fundamental lines {v - 1-0 
centered at /l ~ 4.67 pm) are sensitive to gas at tem perature 
-1500 K. lBlake & Boogertl ( l2004 : lBrittain et al.1 (1200 3) found 
CO gas at r ^ 1000-1500 K in the disk around HD 163296 
and AB Aur. Likewise, CO gas at ~1500 K was detected by 
Ivan den Ancker et all (I2OOII) toward 5 1 Oph. 

The opacity of a dust-poor gas at Tgas = 2000 - 3000 K 
is dominated by H2O and H2, which may eventually deprive 
the inner dust rim of large amount of stellar UV photons below 
1215 A. Indeed, the SED of 51 Oph is well fitted by the emis- 
sion from an optically thick geometrically flat disk and does 
not show significant near-infrared bump, con sistent with the ab- 
sence of a puffed-up rim ( Dullem ond et akfc oOl : Meeus e t alJ 
I2OOII) . 

In summary, the presence of hot CO, traced by the band- 
head emission, and at the same time the absence of inner 
puffed-up rim are probably linked, requiring that r,ot << r^. 
This criterion impl ies simultane ously a high total luminosity 
L - Lt + Lace (.Natta et alJl200lh and a relatively /asf rotating 
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star. Interestingly, this criterion has also been advocated to ex- 
plain t he general absence of i nner puffed-up rim in most T Tauri 
disks jMuzeroUe et alJ2003 ). We compare the thickness of the 
dust free region of 5 1 Oph with two classical Herbig Ae stars 
in Table|2]for grains with radius of 0.1 fim and 1 fim. It is clear 
that 51 Oph constitutes a special case among the Herbig Ae 
stars since it combines a large rotation velocity and a large lu- 
minosity. Indeed CO bandhead emission has not been detected 
toward HD 163296 and AB Aur (R. Waters, private commu- 
nication). The anti-correlation between CO bandhead emission 
and pufFed-up rim should be tested in a larger sample of young 
stars surrounded by a disk. 

Table 2. Comparison of the co-rotation and dust sublima- 
tion radius between 51 Oph and two well studied Herbig Ae 
stars. The spectral type and rotation veloci ty for HD 163296 
and AB Aur are taken from Mora et al .l (l200ll). The other 
stellar chara cteristics are provided by van den Ancker et'aP 
(1199 8'. '2001'). The accretion luminosities were measured by 
iHillen brand et al. (1992). The fits to the Spectral Energy 
Distribution of HD 163296 and AB Aur are discussed in 
lDominiketaD(l2003l) . 





51 Oph 


HD 163296 


AB Aur 


SpT 


B9.5IIIe 


A3Ve 


AOVe+sh 


M. (Me) 


3.8 


2.0+0.2 


2.4+0.2 


L, (Lo) 


260tfo 


26±5 


47+12 


^acc (-^0) 




23 


28 


vsin / km s"^ 


267+5 


133±6 


97+20 


'-.o, (AU) 


0.05 


0.13 


0.33 


ri (AU) (a=1.0/im) 


0.56 


0.18 


0.24 


Td (AU) (a=0.1^m) 


2.50 


0.80 


1.07 


ra - r,„t (AU) (a=1.0/im) 


0.51 


0.05 


-0.09 


rd - not (AU) (a=0.1 /im) 


2.45 


0.38 


0.74 


CO bandhead emission 


yes 


no 


no 


CO fundamental emission 


yes 


yes 


yes 


PufFed-up inner rim 


no 


yes 


yes 



5. Conclusions 

The principle results are summarized here: 

- A high column density of hot molecular (CO) gas is present 
in the inner Astronomical Unit of 5 1 Oph. The gas rotates 
around the star in a Keplerian orbit. 

- The apparent discrepancy on the amount of CO gas around 
51 Oph between the ISO-SWS and FUSE data can be re- 
solved if the disk is geometrically flat and is seen with a 
large angle but not entirely edge-on. 

- The CO bandhead emission is located inside the dust sub- 
limation radius. We propose that only star-disk systems 
where the co-rotation radius is much smaller than the dust 
sublimation radius (r,ot << r^) will show high line over 
continuum contrast CO bandhead emission. This later pre- 
diction can be tested with future high resolving power near- 
infrared interferometric instruments like Amber at the VUT. 



51 Oph is a special case since its rotation speed of 267 km 
s"' is much higher than the average value for Herbig Ae 
stars (80-150 km s"'). 
- Another test for the existence of the dust-free gas would 
be th e detection of strong H? fluorescence lines in the UV 
(e.g.. lHerc7eg et alJll2004 who found H2 UV lines in the 
inner disk around the T Tauri star TW Hya). The obser- 
vation of near-infrared H2 lines with the v=2-l S(l) inten- 
sity at 2.2477 j^m higher than the v= 1-0 S(l) at 2.1218 yum 
would also provide another diagnostic of dust-free gas. 
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